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1 Problem 

In the usual formulation of the Kirchhoff diffraction integral, a scalar field with harmonic 
time dependence at frequency u is deduced at the interior of a charge-free volume from 
knowledge of the field (or its normal derivative) on the bounding surface. In particular, the 
field is propagated forwards in time from the boundary to the desired observation point. 

Construct a time-reversed version of the Kirchhoff integral in which the knowledge of the 
field on the boundary is propagated backwards in time into the interior of the volume. 

Consider the example of an optical focus at the origin for a system with the z axis as the 
optic axis. In the far field beyond the focus a Gaussian beam has cone angle 6*0 = V2ag, and 
the X component of the electric field in a spherical coordinate system is given approximately 
by 

E^{r, e, 0, t) = E(r)e*('=^^--*)e-^'/^o, (1) 

where k = uj/c and c is the speed of light. Deduce the field near the focus. 

Since the Kirchhoff diffraction formalism requires the volume to be charge free, the time- 
reversed technique is not applicable to cases where the source of the field is inside the 
volume. Nonetheless, the reader may find it instructive to attempt to apply the time-reversed 
diffraction integral to the example of an oscillating dipole at the origin. 

2 The Kirchhoff Integral via Green's Theorem 

A standard formulation of Kirchhoff 's diffraction integral for a scalar field with time 
dependence e"""^* is 

^(x) ^ A / ^^(x^Area', (2) 
iTii Js r 

where the spherical waves q^^^'^' -^^^) j f' are outgoing, and r' is the magnitude of vector r' = 
X — x'. 

For a time-reversed formulation in which we retain the time dependence as e"*"^*, the 
spherical waves of interest are the incoming waves ^-'^i^'^'' +^^) j r' . In brief, the desired time- 
reversed diffraction integral is obtained from eq. @) on replacing ihy —i: 

^(x) ^ — / ^4^^(x')rfArea'. (3) 
27r Js r 

For completeness, we review the derivation of eqs. (@)-(0) via Green's theorem. See also, 
sec. 10.5 of ref. M. 
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Green tells us that for any two well-behaved scalar fields and ip, 

J^{<f>^'^^ - ^/'V^0)ciVol = ^(0VV - ^AVV) ■ dS'. (4) 

The surface element dS' is directly outward from surface S. We consider fields with harmonic 
time dependence at frequency u, and assume the factor e"""^*. The wave function of interest, 
ijj, is assumed to have no sources within volume V, and so obeys the Helmholtz wave equation, 

V^V + = 0. (5) 

We choose function ^(x) to correspond to waves associated with a point source at x'. 
That is, 

V2(/) + P(/.= -53(x-x'). (6) 
The well-known solutions to this are the incoming and outgoing spherical waves, 

„±ifcr' 

0±(x,x') = ^, (7) 
where the -|- sign corresponds to the outgoing wave. We recall that 

VV = = -n„, (8) 
where ho points towards the observer at x. Then, 

VV± = Ttkho (l ± ^) (j). (9) 



ikr' 



Inserting eqs. into eq. (H), we find 



1 r 6=^*"'' r / 1 \ 

V'(x) = -— / ■ VV±zA;nJl±— U 

Air Js r' I V tkr / 



dArea', (10) 



where the overall minus sign holds with the convention that n' is the inward normal to the 
surface. 

We only consider cases where the source of the wave ijj is far from the boundary surface, 
so that on the boundary ip is well approximated as a spherical wave. 



^(x')^A , (11^ 



where is the magnitude of the vector = x' — x^ from the effective source point x^ to 
the point x' on the boundary surface. In this case. 



W = ikhs ( 1 ± ^ ) ^, (12) 



where = rg/r^ 
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We also suppose that the observation point is far from the boundary surface, so that 
kr' <^ 1 as well as kvg -C 1. Hence, we neglect the terms in 1/ikr' and l/ikr^ to find 

^(x) = -— / n' ■ (n, ± n,)^(x')c/Area'. (13) 

An Js r' 

The usual formulation, eq. (0), of Kirchhoff 's law is obtained using outgoing waves (+ 
sign), and the paraxial approximation that n' fio ~ hg. The latter tacitly assumes that 
the effective source is outside volume V. 

Here, we are interested in the case where the effective source is inside the volume V, so 
that the paraxial approximation is h' ^ ho ^ — n^. When we use the incoming wave function 
to reconstruct ip{'x, t) from information on the boundary at time t' > t, we use the — sign in 
eq. (|13D to find eq. (1). 

Note that in this derivation, we assumed that ip obeyed eq. (|^) throughout volume V, 
and so the actual source of ip cannot be within V. Our time-reversed Kirchhoff integral (|^) 
can only be applied when any source inside V is virtual. This includes the interesting case of 
a focus of an optical system (sees. 4 and 5). However, we cannot expect eq. (|^) to apply to 
the case of a physical source, such as an oscillating dipole, inside volume V (sec. 6). The laws 
of diffraction do not permit electromagnetic waves to converge into a volume smaller than 
a wavelength cubed, and so eq. (^) cannot be expected to describe the near fields around a 
source smaller than this. 



3 A Plane Wave 

The time-reversed Kirchhoff integral (^ for the x component of the electric field is 

£^2;(obs, now) = — / — — _E'^(r, 6^, 0, future)dArea, (14) 

where r' is the distance from the observation point robs = {x, y, z) in rectangular coordinates 
to a point r = r(sin^cos0, sin^sin0,cos^) on a sphere of radius r in the far field. 
As a first example, consider a plane electromagnetic wave, 

= Eoe'^'''-'^''> = ^Qe*('="™'^^-^*\ (15) 

where the second form holds in a spherical coordinate system (r, 6, 0) where 6 is measured 
with respect to the z axis. We take the point of observation to be {x,y,z) = (0,0, ro), and 
evaluate the diffraction integral (p^ over a sphere of radius r ^ tq. In the exponential 
factor in the Kirchhoff integral, we approximate r' as 

r' ^ r — r ■ Fobs = r — Tq cos 9, (16) 

while in the denominator we approximate r' as r. Then, 

•L ^1 ^27r p-ik{r-ro cos 9) 

E^(obs) ^ — / r^d COS 9 dcj) Eod^'^"'^ 

2tt J-1 Jo r 



r + ro 
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where we ignore the rapidly oscillating term e ^''('^^'+^0) g^g unphysical. 

This verifies that the time- reversed diffraction formula works for a simple example. 



4 The Transverse Field near a Laser Focus 

We now consider the far field of a laser beam whose optic axis is the z axis with focal 
point at the origin. The polarization is along the x axis, and the electric field has Gaussian 
dependence on polar angle with characteristic angle ^0 ^ 1- Then, we can write 

E,{r, e, 0) = E(r)e*'-'-e-^'/^o, (18) 

where E{r) is the magnitude of the electric field on the optic axis at distance r from the 
focus. In the exponential factor in the Kirchhoff integral (|1^, r' is the distance from the 
observation point r\rmobs = [x, y, z) to a point r = r(sin 9 cos 0, sin 9 sin 0, cos 9) on the sphere. 
We approximate r' as 

r' r — i ■ Fobs = r — X sin 9 cos (p — y sin 9 sin cf) — z cos 6', (19) 

while in the denominator we approximate r' as r. Inserting eqs. ([T8|) and (|1^) into (p^), we 
find 

^^(obs) = ^^^-^(^) 1^ g»fc^coseg-eveg^cosg Z'^" sine cos <A+»fc^ sine sin 
27r J-i Jo 

= ikrE{r) y\*'="^°^^e-^'/^oJo(A;psin^)rfcos^, (20) 

where 



p=^x^ + y\ (21) 

and Jo is the Bessel function of order zero. 

Since we assume that the characteristic angle ^0 of the laser beam is small, we can 
approximate cos 6' as 1 — 6*^/2 and kpsin9 as kp9. Then, we have 



E,(obs) ^ zA;rE(r)e''="y e-(2/^«+''^")^'/Vo (^V2A;py^^rf(eV2) 

2{1 + ik9iz/2) ' 



(22) 



where the Laplace transform, which is given explicitly in [^, can be evaluated using the 
series expansion for the Bessel function. This expression can be put in a more familiar form 
by introducing the Rayleigh range (depth of focus). 



2 

k9i 

and the so-called waist of the laser beam, 



^0 = 117^2' (23) 



2 

Wo = 9oZo = — . (24) 

KtJQ 
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We define the electric field strength at the focus {p = 0, z = 0) to be Eq, so we learn that 
the far-field strength is related by 

E{r) = -z^Eo. (25) 



The factor —i = e ^^"^ is the 90° Guoy phase shift between the focus and the far field. Then, 
the transverse component of the electric field near the focus is 




Figure 1: The electric field E^{x, 0, z) of a linearly polarized Gaussian beam 
with diffraction angle 9o = 0.45. 

The Gaussian beam (^6|) could also be deduced by a similar argument using eq. (|^), 
starting from the far field of the laser before the focus. The form (^) is symmetric in z 
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except for a phase factor, and so is a solution to the problem of transporting a wave from 
z = —r to z = +r such that the functional dependence on p and z is invariant up to a phase 
factor. One of the earliest derivations ^ of the Gaussian beam was based on the formulation 
of this problem as an integral equation for the eigenfunction (|26D . 



5 The Longitudinal Field 

Far from the focus, the electric field E(r) is perpendicular to the radius vector r. For a field 
linearly polarized in the x direction, there must also be a longitudinal component related 
by 

E-r = E^.sin^cos0 + E^cos^ = O. (27) 

Thus, far from the focus, 

E, (r) = -E,^ (r) tan 6 cos (f). (28) 
Then, similarly to eqs. (|l^) and (pO[), we have 



ikr' 



E,(obs) = ^ f ^—E,{r)dATea 
Ztt J r 

^ _ lkrE{r) /■! gifc.cos0g-0Vegtan^rfcOS^ sin e cos 0+^., sin e sin <^ ^ 

27r J-1 Jo 

^ _ tkxzoEo e»fc.cos0g-eVegtangJi(A;psing)dcosg, (29) 
p J-i 

using eq. (3.937.2) of §. 

We again note that the integrand is significant only for small 6, so we can approximate 
eq. (p9D as the Laplace transform 



E,{x,y,z) ^ -texzoEoe''''V2 e-^^/^o+^'^^^'^'/^y^j^ J^^^^^^j ^(^2/3) 
ik'^9^xzoEoe'^'e-p'/'"o(^+''/'°^ 



A{l + iz/zo) 
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X E^{x,y,z) 

= -iOq— , ■ / 30) 

Wo {1 + iz/zo) 

with Ex given by eq. (pG)) . Figure ^ plots this field. 

Together, the electric field components given by eqs. (^) and ( POD satisfy the Maxwell 
equation V ■ E = to order 6^ i, 0, 



6 Oscillating Dipole at the Origin 

We cannot expect the Kirchhoff diffraction integral to apply to the example of an oscillating 
dipole, if our bounding surface surrounds the dipole. Let us see what happens if we try to 
use eq. (H) anyway. 
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Figure 2: The electric field Ez{x, 0, z) of a linearly polarized Gaussian beam 
with diffraction angle 6q = 0.45. 



The dipole is taken to be at the origin, with moment p along the x axis. Then, the x 
component of the radiation field is 

Ex = k'^p sinOx . (31) 

r 

where 6x is the angle between the x axis and a radius vector to the observer. We consider 
an observer near the origin at (x, y, z) = (0, 0, tq), for which sin^a; = 1, and so 

E^(obs) = k^p^^. (32) 

We now attempt to reconstruct this field near the origin from its value on a sphere of 
radius r using the time- reversed Kirchhoff integral (|^). We use a spherical coordinate system 
(r, 6, 0) that favors the z axis. Then, the x component of the radiation field on the sphere 
of radius r is ^ 

Ex{r, 9, (p) = k'^p^l-sin^Ocos^cp^. (33) 

This form cannot be integrated analytically, so we use a Taylor expansion of the square root, 
which will lead to an expansion in powers of l/rg. It turns out that the coefficient of the 
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1/ro term, which is our main interest, is very close to that if we simply approximate the 
square root by unity. For brevity, we write 

E,{r,e,(j))^k^p—. (34) 

In the time- reversed Kirchhoff integral (^, we make the usual approximation that r' = 
r — TqCOS^ in the exponential factor, but r' = r in the denominator. Then, using eq. (0) 
we have 

? Jc^ TIP ^^'^ /"l /'27r pihv 

EJohs) ^ / r^dcosO / d0e^'='-«'=°^^ — 

27rr J-i Jo r 

gjfcro g-ifcro 

= k'^p k'^p 

, n sin krn , , 

= lik^p— -. 35 

Zero 

The first, outgoing wave in middle line of eq. (|35[) is the desired form, but the second, 
incoming wave is of the same magnitude. Together, they lead to the form •sm.ik.ro) jkr^ 
which is nearly constant for kr^ < 1. The presence of outgoing as well as incoming waves is 
to be expected because dipole radiation is azimuthally symmetric about the x axis. In the 
absence of a charged source at the origin, an outgoing wave at ^ = vr must correspond to an 
incoming wave at ^ = 0. 

The result that the reconstructed field is uniform for distances within a wavelength of the 
origin is consistent with the laws of diffraction that electromagnetic waves cannot be focused 
to a region smaller than a wavelength. Far fields of the form (^TJ) could only be propagated 
back to the form of dipole fields near the origin with the addition of nonradiation fields tied 
to a charge at the origin. Such a construction is outside the scope of optics and diffraction. 
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